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ABSTRACT 


This  annual  report  details  progress  in  basic  research  in  thermoacoustic  heat 
transport  made  during  the  period  October  1,  1992  through  September  30,  1993.  Our 
major  research  efforts  in  FY  1993  were  concentrated  in  five  areas:  1)  an  analysis  of  the 
initial  buildup  of  oscillations  in  thermoacoustic  prime  movers;  2)  an  investigation  of 
stability  curves  for  a  thermoacoustic  prime  mover,  3)  initial  measurements  of  the  velocity 
field  in  thermoacoustic  engines  using  laser  Doppler  anemometry,  4)  measurements  of 
heat  exchanger  performance  in  a  thermoacoustic  prime  mover,  and  5)  a  preliminary 
investigation  of  heat  driven  refrigerators.  In  addition  to  these  areas,  FY  1993  research 
efforts  also  included  a  collaborative  experiment  on  traveling  wave  thermoacoustic  effects. 
A  brief  review  of  the  motivation  behind  these  studies  and  a  summary  of  the  important 
results  are  contained  in  this  technical  report.  A  publications,  patents,  presentations,  and 
honors  report  is  also  included. 


BACKGROUND 

The  central  theme  of  this  research  project  has  been  that  thermoacoustic  heat 
transport  phenomena  are  not  well  understood  at  high  acoustic  amplitudes.  That  this 
central  theme  is  still  true  has  been  most  recently  confirmed  in  Ref.  1.  In  that  work,  Swift 
presents  an  extensive  set  of  measurements  of  the  performance  of  a  large  thermoacoustic 
prime  mover.  He  uses  a  steady  state,  energy  balance  technique  based  on  linear  theory  to 
calculate  the  steady  state  parameters  of  his  engine.  A  heat  input  is  specified,  in  watts,  for 
the  system.  This  heat  input  is  converted  to  an  acoustic  power  output  according  to 
thermoacoustic  theory.  Specifying  the  input  in  real  units  yields  a  dimensional,  rather  than 
nondimensional,  calculated  acoustic  output.  Swift  observes  significant  discrepancies 
between  the  predicted  and  measured  performance,  especially  at  high  amplitudes. 
Pippard2  points  out  that  the  factors  that  determine  the  steady  state  amplitude  of  a 
spontaneously  oscillating,  such  as  a  prime  mover,  may  ultimately  be  nonlinear. 
Therefore,  it  is  not  surprising  that  discrepancies  are  observed.  These  nonlinearities  may 
be  attributed  to  different  physical  phenomena  associated  with  prime  movers.  Three 
candidates  from  a  nonexhaustive  list  immediately  come  to  mind.  They  are  harmonic 
generation,  velocity  dependent  effects,  and  heat  exchanger  performance. 

We^  and,  independently,  Swift1  have  addressed  harmonic  generation  previously. 
The  current  level  of  understanding  of  the  role  of  harmonic  generation  in  high  amplitude 
acoustic  standing  waves  is  generally  sufficient  for  us  to  mitigate  its  effects  in 
thermoacoustics.  Therefore,  we  turn  to  other  limiting  factors.  Experimental  evidence 
that  some  acoustic  velocity  dependent  effect  (such  as  acoustic  streaming,  vortex 
shedding,  or  turbulence)  may  influence  acoustic  heat  transport  was  pointed  out  early  on  in 
our  research.4  To  investigate  velocity  dependent  effects  further,  we  have  initiated  laser 
Doppler  anemometer  studies  of  the  velocity  field  in  thermoacoustic  engines. 


The  role  of  heat  exchange  is  paramount  in  thermoacoustics  and  there  are  a  number 
of  basic  research  issues  related  to  it.  There  is  evidence  that  the  limiting  factor  accounting 
for  the  observed  steady  state  amplitudes  in  prime  movers  is  inadequate  heat  exchanger 
performance.5  One  manifestation  of  inadequate  heat  exchangers  is  the  finite  capacity  for 
transferring  heat  into  the  engine  from  an  external  heat  source.  Limitations  in  available 
power  input  may  not  seem  at  first  to  be  a  nonlinear  mechanism.  However,  it  is  in  the 
same  sense  that  a  finite  power  supply  output  voltage  level  provides  the  nonlinearity  that 
limits  the  amplitude  of  a  spontaneously  oscillating  electrical  circuit.  The  growth  rate  of 
the  transient  buildup  of  oscillations  changes  abruptly  when  the  circuit  output  voltage 
reaches  the  power  supply  limits  -  the  output  clips.  As  long  as  the  thermal-to-acoustic 
conversion  involves  only  linear  processes,  the  linear,  steady  state  energy  balance 
approach  discussed  earlier  works.  However,  if  nonlinear  processes  are  involved,  the 
calculations  may  have  to  be  altered  to  take  into  account  the  transient  part  of  the  solution. 
Knowing  at  what  acoustic  levels  various  processes  become  important  then  becomes 
critical  knowledge. 

If  it  turns  out  that  heat  exchanger  performance  is  the  limiting  factor,  then  only 
after  the  limitations  are  eliminated  can  other,  perhaps  more  fundamental,  limiting  factors 
be  investigated.  If  heat  exchanger  performance  is  not  now  the  prime  limiting  factor,  then 
it  is  likely  to  become  so  when  the  current  limiting  factors  are  better  understood  and  their 
effects  lessened  or  when  required  heat  loads  increase.  Currently,  thermoacoustic 
refrigerators  handle  heat  loads  of  the  order  102  W.  However,  heat  loads  of  the  order  104 
W  are  required  for  Navy  refrigeration  and  cooling  purposes.  Heat  loads  of  this 
magnitude  place  thermoacoustics  in  a  completely  unexplored  regime.  We  are  currently 
investigating  heat  exchanger  performance  in  high  amplitude  sound  fields  generated  by  a 
prime  mover.  This  work  in  performed  in  collaboration  with  Prof.  Thomas  Hofler.  It  is 
expected  that  knowledge  gained  from  prime  mover  experiments  will  be  directly 
transferable  to  refrigerator  technology. 


This  background  has  centered  around  prime  movers.  However,  it  is  important  to 
keep  in  mind  that  much  of  what  is  learned  about  prime  movers  can  be  carried  over  to 
improve  thermoacoustic  refrigeration  and  cooling  technology. 


SUMMARY  OF  PROGRESS 


1)     Analysis  of  the  Initial  Buildup  of  Oscillations  in  a  Thermoacoustic  Prime  Mover 

As  discussed  earlier,  predicting  the  steady  state  amplitude  of  acoustic  oscillations 
in  prime  movers  is  one  of  the  fundamental  questions  yet  to  be  fully  answered  about 
thermoacoustics.  The  purpose  of  this  research  is  to  address  some  issues  related  to  the 
transition  to  steady  state  above  onset.  In  particular,  we  analyze  the  initial  buildup  of 
oscillations  in  a  prime  mover  and  track  the  quality  factor  Q  of  a  prime  mover  through 
onset.  Full  details  of  this  work  can  be  found  in  Ref.  6. 

The  transition  to  onset  in  a  prime  mover  is  conveniently  discussed  in  terms  of  its 
Q.  The  acoustic  power  W  generated  by  a  prime  mover  can  be  related  to  Q  through  the 
relationship 

Q  =  -  —f1.  (l) 

W 

A  typical  prime  mover  is  comprised  of  five  sections:  the  ambient  duct,  the  ambient  heat 
exchanger,  the  prime  mover  stack,  the  hot  heat  exchanger,  and  the  hot  duct.  W  is  the 
sum  of  the  power  output  of  the  five  individual  sections  of  the  prime  mover, 

W  =  W     .  +  W    .  .    +  W      +  W.  , .    +  W.  ,.  (2) 

amb         ambhx         pms         hot  hx  hot  v  ' 


The  subscripts  amb,  amb  hx,  pms,  hot  hx,  and  hot  refer  to  the  ambient  duct,  the  ambient 

heat  exchanger,  the  prime  mover  stack,  the  hot  heat  exchanger,  and  the  hot  duct, 

respectively.     W   is  defined  such  that  it  is  positive  when  power  is  generated.    Q  is 

typically  expressed  in  terms  of  dissipated  power,  which  is  defined  so  as  to  be  positive 

when  power  is  lost.  Because  we  are  interested  in  power  generation,  we  use  W  and 
introduce  a  minus  sign  in  Eq.  (1).  Kj  is  the  stored  acoustic  energy  and  co  is  the  angular 

frequency  of  the  oscillations. 

All  the  terms  on  the  right-hand  side  of  Eq.  (2)  are  inherently  negative 
(representing  dissipation)  except  for  W      .   W       can  be  either  positive  or  negative, 

depending  on  the  temperature  difference  AT  imposed  across  the  stack  and  its  position  in 
the  standing  wave.  As  AT  increases  from  zero,  the  thermal  losses  in  the  stack  decrease. 
Hence,  the  net  losses  in  the  prime  mover,  consisting  of  thermal  and  viscous  wall  losses  in 
each  of  the  five  sections  (neglecting  attenuation  in  the  gas)  decrease.  The  Q  increases. 
As  AT  increases  further,  the  thermal  losses  in  the  stack  decrease  to  zero  and  then  become 
negative,  representing  gain.  At  some  higher  AT,  the  thermal  gain  in  the  stack  exactly 
balances  the  losses  in  the  stack  and  W  =  0.  However,  W  is  still  negative,  due  to  the 
other  four  terms  in  Eq.  (2),  and  so  Q  is  positive  and  finite.  At  a  sufficiently  high  AT,  the 
thermal  gain  overcomes  all  the  losses  andlV  =  0.  At  this  point  Q  is  infinite,  the  net  loss 
in  the  prime  mover  zero.  The  prime  mover  is  in  a  state  of  unstable  equilibrium.  Any 
infinitesimal  increase  in  AT  will  make  W       >  0,  the  net  loss  and  Q  negative;  the  prime 

mover  will  go  into  spontaneous  self-oscillation.  The  complete  evolution  of  Q  as  just 
described  has  never  been  previously  demonstrated  in  the  literature  for  an  individual  prime 
mover. 

Regarding  the  question  of  what  determines  the  steady  state  amplitude  above 
onset,  Pippard  states  that  "non-linearity  is  ...  an  ultimately  inevitable  feature  of  a 
spontaneously  oscillating  system."2  The  implication  is  that  the  factors  determining  the 


steady  state  amplitude  in  a  prime  mover  may  not  be  contained  in  linear  thermoacoustic 
theory.  If  one  uses  linear  acoustic  theory  to  derive  an  expression  for  Q,  one  finds  that  the 
result  is  independent  of  the  acoustic  pressure  amplitude.  Therefore,  according  to  linear 
theory,  once  Q  is  finite-negative  it  will  remain  so  regardless  of  the  pressure  amplitude  in 
the  prime  mover.  Therefore,  because  Q  =  «»  at  steady  state,  steady  state  can  never  be 
reached  according  to  linear  theory.  However,  steady  state  oscillations  are  observed  in 
prime  movers.  Therefore,  the  states  (Q  =  ^Oonset  and  (Q  =  Steady  state  are  fundamentally 
different. 

One  step  along  the  path  to  understanding  the  factors  that  govern  steady  state  is  to 
analyze  the  initial  buildup  of  oscillations  in  a  prime  mover  above  onset.  During  the 
initial  buildup,  the  acoustic  amplitudes  are  low  and  most  of  the  complicating  factors  can 
be  neglected.  Analysis  of  the  buildup  and,  in  particular,  at  what  amplitudes  deviations 
between  theory  and  measurement  begin  may  help  to  identify  particular  mechanisms. 

Previously,7  we  have  discussed  the  evolution  of  Q  up  to  onset.  This  work  extends 
our  investigation  of  prime  movers  through  onset.  The  present  measurements,  coupled 
with  previous  ones  (reported  in  Ref.  7)  show  the  complete  evolution  of  Q  .  We  also  find 
that  the  initial  buildup  of  oscillations  above  onset  is  consistent  with  linear  theory.  It  is 
shown  that  the  results  are  consistent  with  predictions  based  on  a  standing  wave  analysis 
of  thermoacoustics.7 


2)    Stability  Curves  for  a  Thermoacoustic  Prime  Mover 

The  onset  of  self-oscillation  in  prime  movers  can  be  considered  in  terms  of  a 
transition  from  a  stable  state  to  an  unstable  state.  The  stable  state,  below  onset,  is  one  for 
which  the  prime  mover  as  a  whole  has  a  net  positive  damping.  Any  oscillations  excited 
in  a  prime  mover  below  onset  decay  in  time.  The  unstable  state,  above  onset,  is  one  for 


which  the  damping  is  negative.  Oscillations  at  the  resonance  frequency  initially  exhibit 
exponential  growth  in  amplitude.  Eventually  steady  state  is  obtained  and  the  growth  rate 
is  zero.  The  steady  state  acoustic  pressure  amplitude  is  typically  on  the  order  of  a  few 
percent  or  more  of  the  mean  gas  pressure. 

For  a  given  prime  mover  geometry,  mode  and  gas  type,  the  temperature  difference 
AT  required  for  onset  ATonset  is  a  function  of  the  mean  gas  pressure  Pm.  A  plot  of  ATonset 
versus  Pm  defines  the  so  called  stability  curve  for  a  prime  mover.  Two  such  curves  are 
shown  in  Fig.  1.  The  dashed  (dash-dot)  curve  corresponds  to  the  fundamental  (second) 
mode.  In  order  for  self-sustained  oscillations  to  exist,  both  parameters  ATonset  and  Pm 
must  lie  within  the  region  of  instability,  above  the  curves  in  Fig.  1. 

An  extensive  series  of  investigations  on  the  behavior  of  simple  prime  movers 
(although  they  didn't  call  them  that)  was  conducted  by  Yazaki,  et  al.8-10  Their  prime 
movers  consisted  of  U-shaped  stainless  steel  tubes  (inner  radius  1.2  -  4.7  mm,  length  1.5 
m  and  greater).  The  temperature  of  the  lower  part  of  the  U  was  adjusted  by  a  continuous 
flow  of  cold  helium  gas.  The  upper  p^rt  of  the  U  was  maintained  at  room  temperature 
(293  K).  They  investigated  both  the  limits  of  stability  of  their  system  as  well  as  the 
highly  nonlinear  nature  of  the  oscillations  above  onset.  They  discovered  a  wealth  of 
complex  behavior,  including  phase  locking,  mode  competition,  quasiperiodicity  and 
chaos. 

We  have  undertaken  a  similar  investigation  on  a  practical  prime  mover.  In 
particular,  we  have  concentrated  on  the  excitation  of  the  second  longitudinal  mode  of  the 
prime  mover  and  the  simultaneous  excitation  of  both  the  fundamental  and  second  modes. 
The  region  of  two  mode  excitation  in  practical  prime  mover  has  received  little  attention 
previously.   The  results  are  summarized  below.   Full  details  can  be  found  in  Ref.  11. 
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The  measured  stability  curves  for  the  fundamental  and  second  modes  of  the  prime 
mover  are  indicated  by  the  symbols  in  Fig.  1.  The  o's  indicate  the  measured  onset 
conditions  for  the  fundamental  mode.  The  *'s  indicate  the  measured  onset  conditions  for 
the  second  mode.  The  +'s  indicate  the  onset  conditions  for  the  fundamental  mode,  after 
the  second  mode  is  above  onset.  These  last  data  points  were  obtained  by  fixing  AT  and 
increasing  Pm.  We  estimate  that  the  maximum  error  is  ±5  K  for  ATonset  and  ±1  to  2%  for 
Pm.  It  is  seen  that  the  predicted  and  measured  stability  limits  are  in  reasonable  agreement 
for  the  fundamental  mode  at  most  mean  pressures.  The  predicted  trend  for  the  second 
mode  is  consistent  with  the  measurements,  yet  the  predicted  values  of  AT  onset  are  higher 
than  those  measured.  The  reason  for  this  discrepancy  with  the  second  mode  is  not 
known,  but  is  consistent  with  results  reported  in  Ref.  2. 

There  are,  however,  two  more  interesting  aspects  of  Fig.  1.  The  first  is  that  the 
slope  of  the  measured  stability  curve  for  the  fundamental  is  positive  when  the  second 
mode  is  excited.  (These  points  are  represented  by  the  +'s.)  The  predicted  slope  of  the 
stability  curve  at  low  pressures  is  negative.  It  is  not  completely  surprising  that  the  slope 
is  different  than  that  predicted  because  the  theoretical  curves  do  not  take  into  account  the 
presence  of  other  modes.  Apparently  the  presence  of  the  second  mode  does,  in  fact, 
affect  the  stability  of  the  fundamental.  Note  that  the  slope  of  the  measured  low  pressure 
stability  limit  of  the  second  mode  is  negative,  as  predicted.  Also  note  that  the  second 
mode  is  the  only  one  present  in  this  region.  In  other  words,  there  is  no  other  mode 
present  to  interfere  with  it.  One  possible  interaction  is  that  the  presence  of  the  second 
mode  alters  the  temperature  distribution  in  the  stack  significantly,  so  that  AT  is 
unrepresentative  of  the  actual  conditions  in  the  stack.  Another  possible  explanation  is 
that  there  is  nonlinear  mode  competition  between  the  fundamental  and  second  mode. 
When  the  second  mode  is  present,  it  has  the  advantage  over  the  fundamental  well  into  the 
region  where  the  fundamental  should  be  unstable.  The  behavior  of  the  stability  curve  of 


the  fundamental  at  low  pressures  shows  that  below  approximately  140  kPa,  the 
fundamental  mode  will  never  be  excited,  regardless  of  how  high  AT  becomes!  This  result 

is  unexpected. 

It  is  also  observed  that  one  mode  can  suppress  the  other.  This  is  the  second 
interesting  aspect  of  Fig.  1  mentioned  above.  The  measured  stability  curve  for  the  second 
mode  does  not  extend  to  the  right  of  the  intersection  with  that  of  the  fundamental.  In 
other  words,  we  never  observe  onset  of  the  second  mode  once  the  fundamental  is  excited. 
Two  methods  were  used  to  search  for  second  mode  onset.  In  one,  Pm  was  fixed  and  AT 
increased  up  to  a  maximum  value  of  approximately  550  K,  well  above  that  which  should 
be  required  for  second  mode  onset.  In  the  other  case,  we  fixed  AT  and  increased  Pm.  In 
neither  case  did  the  second  mode  reach  onset. 

Equally  interesting  is  what  happens  to  the  second  mode  when  the  fundamental 
mode  reaches  onset.  To  use  a  specific  example,  if  AT  is  set  to  450  K  and  Prn  is  increased 
from  zero,  the  second  mode  reaches  onset  at  a  pressure  of  approximately  115  kPa.  As  Pm 
is  increased  further,  the  amplitude  of  the  second  mode  increases.  When  the  mean 
pressure  reaches  approximately  150  kPa,  we  see  the  fundamental  mode  begin  to  grow. 
For  a  brief  period,  both  modes  are  simultaneously  excited  with  the  amplitude  of  the 
fundamental  increasing  and  that  of  the  second  mode  decreasing.  Soon,  the  second  mode 
is  completely  gone  and  only  the  fundamental  remains.  Further  increase  of  Pm  only 
results  in  the  growth  of  the  fundamental.  The  second  mode  never  comes  back,  even 
though  it  should  be  unstable. 

To  investigate  this  behavior  further,  we  decided  to  selectively  inhibit  the 
fundamental  mode,  without  affecting  the  second.  This  is  accomplished  by  inserting  an 
annular  disk  in  the  prime  mover  at  the  location  of  the  velocity  antinode  of  the 
fundamental  mode.  The  3.13-mm  thick  brass  disk  has  a  18.3  mm  diameter  hole  in  the 
center.  The  idea  is  that  the  reduced  cross  section  presents  a  significant  flow  impedance  to 
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the  fundamental  mode,  thereby  raising  its  ATonset.  But,  because  the  disk  is  near  the 
velocity  node  of  the  second  mode,  the  second  mode  onset  conditions  will  be  affected  very 
little.  By  inhibiting  the  fundamental,  the  fundamental's  influence  on  the  second  mode 
should  be  diminished. 

The  stability  curves  for  the  modified  prime  mover  are  shown  in  Fig  2.  The 
symbols  have  the  same  meaning  as  in  Fig.  1.  The  solid  lines  simply  connect  the  data 
points  and  serve  as  guides  to  the  eye.  Several  points  are  immediately  evident.  First,  the 
left  branch  of  the  stability  curve  for  the  second  mode  is  unaltered  by  the  presence  of  the 
disk.  However,  the  stability  curve  for  the  fundamental  shifts  toward  higher  AT  and  Pm. 
In  addition  to  shifting  the  stability  curve,  the  disk  also  results  in  a  much  lower 
fundamental  amplitude.  This  apparently  allows  the  second  mode  to  dominate  over  the 
fundamental  to  a  larger  extent,  as  indicated  by  region  II  occupying  more  parameter  space 
in  Fig.  2  than  in  Fig.  1.  Second,  the  stability  curve  for  the  second  mode  now  extends  into 
the  region  where  the  fundamental  mode  previously  dominated.  The  reduced  fundamental 
amplitude  apparently  prohibits  it  from  quenching  the  second  mode.  In  the  region  labeled 
in,  both  modes  are  simultaneously  excited  and  coexist  in  steady  state. 

We  recorded  waveforms  for  the  modified  prime  mover  at  AT  equal  to  450  K  and 
mean  gas  pressure  of  150,  240  and  440  kPa.  These  parameters  place  us  in  regions  I,  II 
and  III  of  Fig.  2.  The  waveforms  for  the  points  (Pm,  AT)  =  (150  kPa,  450  K)  and  (440 
kPa,  450  K)  show  that  only  one  mode  (the  second  and  fundamental,  respectively)  is 
excited.  The  situation  is  different  at  the  point  (240  kPa,  450  K)  in  region  III.  In  this 
region  both  modes  are  excited  simultaneously.  Spectral  analysis  reveals  that  the 
frequency  of  the  fundamental  mode  is  484.0  Hz  and  that  of  the  second  mode  is  996. 1  Hz. 
Because  the  fundamental  and  second  modes  are  not  harmonic,  the  second  mode  is  easily 
distinguished  from  harmonics  of  the  fundamental  mode.  More  importantly,  it  means  that 
the  combination  is  not  periodic.  It  is  quasiperiodic. 
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The  results  of  this  work  point  to  several  areas  for  further  work.  First,  there  is  need 
for  a  model  that  accurately  predicts  the  stability  curves  by  accounting  for  the  interaction 
between  modes.  Also,  measurements  of  the  stability  curves  in  parameter  regions  where 
both  modes  exist  may  prove  to  be  a  useful  testbed  for  various  thermoacoustics  models. 
The  use  of  a  series  of  discs  with  progressively  larger  holes  and  progressively  larger 
fundamental  amplitudes  might  be  useful  in  studying  the  reported  effects  in  more  detail. 
Finally,  a  more  extensive  investigation  of  (Pm,  AT)  parameter  space  may  reveal  regions 
of  more  complex  behavior.  Such  a  study  may  lead  to  novel  applications  regarding  the 
control  of  oscillations  in  prime  movers. 


3)    Initial  Measurements  of  the  Velocity  Field  in  Thermoacoustic  Engines  Using 
Laser  Doppler  Anemometry 

As  discussed  in  the  Background  section,  we  have  initiated  an  investigation  of  the 
velocity  field  in  thermoacoustic  engines  using  a  laser  Doppler  anemometer  (LDA). 
These  preliminary  measurements  consist  of  probing  the  velocity  field  in  an  acoustic 
resonator  driven  with  a  high  intensity  compression  driver.12'13  The  first  set  of 
measurements  consisted  of  measuring  the  acoustic  velocity  amplitude  at  a  velocity 
antinode  as  a  function  of  the  acoustic  pressure  amplitude.  The  velocities  were  converted 
to  pressures  using  the  standard  value  of  specific  acoustic  impedance  for  air  (415  kg  m*2  s" 
1).  A  graph  of  the  so  called  "derived  pressure"  versus  the  actual  measured  pressure  is 
shown  in  Fig.  3.  The  solid  line  is  a  least  squares  fit  to  the  data.  The  slope  of  the  fit  is 
0.997.  Data  was  obtained  for  pressures  ranging  from  20  to  4000  Pa.  (4000  Pa  was  about 
the  most  we  could  get  from  our  driver  without  blowing  the  voice  coil.)  We  also 
measured  the  velocity  as  a  function  of  position  along  the  axis  of  the  resonator  at  a  fixed 
pressure    amplitude.      The    results    show    the    expected    sinusoidal    variation. 
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Figures  4  and  5  show  the  acoustic  velocity  and  acoustic  pressure  as  functions  of 
time  for  two  pressure  amplitudes.  It  is  seen  that  at  low  amplitudes  the  two  waveforms 
overlap.  This  is  expected  if  one  is  in  the  linear  acoustic  regime.  At  higher  amplitude, 
Fig.  5,  harmonics  appear  in  both  waveforms.  Analysis  of  the  amplitudes  of  the 
harmonics  should  yield  the  nonlinear  coefficient  of  the  gas. 

Finally,  we  put  a  simple  stack  constructed  from  microscope  slides  in  the  resonator 
and  made  measurements  at  various  distances  from  the  stack  entrance.  The  results  are 
shown  in  Fig.  6  for  three  different  pressure  amplitudes.  Negative  positions  are  outside 
the  stack,  the  edge  is  at  zero  and  positive  positions  are  inside  the  stack.  The  dashed  lines 
represent  the  change  in  velocity  expected  from  the  change  in  cross  section  and  position  in 
the  standing  wave. 

These  preliminary  measurements  indicate  that  LDA  is  a  useful  tool  in  probing  the 
velocity  field.  They  also  show  that,  at  least  in  the  parameter  range  covered,  no  significant 
perturbations  in  fluid  velocity  on  acoustic  time  scales  are  apparent.  We  have  made  a 
video  tape  of  flow  processes  in  the  resonator  which  shows  significant  activity,  albeit  at 
slower  time  scales.  We  intend  to  investigate  these  time  scales  more  in  the  future. 


4)     Measurements  of  Heat  Exchanger  Performance  in  a  Thermoacoustic  Prime 
Mover 

This  research  is  being  done  in  collaboration  with  Prof.  Thomas  Hofler  and  Dr. 
David  Gardner.  To  assess  heat  exchanger  performance  we  have  built  a  modular  prime 
mover  which  allows  installation  of  different  heat  exchangers  without  altering  the 
geometry  of  the  remainder  of  the  prime  mover.  The  prime  mover  is  designed  to  operate 
between  room  temperature  and  liquid  nitrogen  temperatures.  The  working  substance  is 
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neon.  The  hot  end  of  the  prime  mover  is  fitted  with  a  movable  piston  which  allows  us  to 
vary  the  position  of  the  thermoacoustic  elements  (stack  and  heat  exchangers)  relative  to 
the  standing  wave.  The  prime  mover  is  instrumented  with  several  thermocouples  which 
permit  measurement  of  the  temperatures  of  the  hot  and  cold  ends  of  the  resonator,  the 
temperature  of  the  hot  and  cold  heat  exchangers  both  near  the  wall  of  the  resonator  and 
near  the  center  of  the  heat  exchanger  (along  the  axis  of  the  resonator),  and  the 
temperatures  of  the  hot  and  cold  ends  of  the  stack  near  the  stack/heat  exchanger  junction. 
Knowledge  of  the  temperature  distribution  in  the  transverse  direction  in  the  heat 
exchanger  (perpendicular  to  the  resonator  axis)  allows  calculation  of  the  heat  flow  into  or 
out  of  the  heat  exchanger.  Knowledge  of  the  temperature  difference  across  the  stack/heat 
exchanger  junction  provides  information  about  the  acoustic  heat  transport  across  the 
junction. 

We  have  completed  an  initial  set  of  measurements.14'1^  The  important  results 
are:  1)  acoustic  pressure  amplitudes  as  high  as  29%  of  mean  gas  pressure  were  generated 
in  the  prime  mover;  2)  short  heat  exchangers  performed  better  than  might  have  expected 
from  simple,  boundary  layer  arguments;  3)  peak-to-peak  displacement  amplitudes  as 
much  as  10  times  the  heat  exchanger  length  were  achieved;  and  4)  peak-to-peak 
displacement  amplitudes  equal  to  the  stack  length  (2.5  cm)  were  achieved.  Such  large 
displacements  are  not  easily  handled  by  standard  thermoacoustic  theory,  which  assumes 
small  displacements  compared  to  the  stack  length. 


5)     Preliminary  Investigation  of  Heat  Driven  Refrigerators 

Investigation  of  heat  driven  thermoacoustic  refrigerators  has  been  one  of  the  long 
term  goals  of  our  research  program.  The  development  of  a  refrigerator  is  certainly  more 
applied  than  most  of  our  previous  work.  But,  because  a  heat  driven  refrigerator  consists 
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of  both  a  prime  mover  and  a  heat  pump,  it  constitutes  a  much  more  demanding  test  of  our 
understanding  of  thermoacoustics.  Although,  Wheadey  built  a  prototype  device  (with 
which  he  was  able  to  achieve  cold  end  temperatures  below  0  °C)  shortly  before  his  death, 
there  has  been  little  subsequent  work  on  this  type  of  device. 

During  FY  1993,  we  began  a  collaboration,  through  an  IPA,  with  Dr.  Pat  Arnott  at 
the  Desert  Research  Institute  to  adapt  his  thermoacoustic  modeling  techniques  to  handle 
heat  driven  refrigerators.  At  this  point,  Dr.  Arnott  has  finished  writing  the  code  and  has 
done  some  optimization  analysis16 


6)     Traveling  Wave  Thermoacoustic  Effects 

Colleagues  at  the  University  of  Mississippi  (Profs.  Richard  Raspet  and  Henry 
Bass)  have  become  interested  in  traveling  wave  thermoacoustic  effects.  They  developed 
a  theory  that  shows  that  the  traveling  wave  phasing  between  pressure  and  velocity  leads 
to  heat  transport  terms  not  contained  in  standing  wave  thermoacoustic  theory.17 
Previously,  we  have  considered  the  possibility  of  using  a  thermoacoustic  device  to 
terminate  a  traveling  wave  tube.18  This  mutual  interest  in  traveling  wave  effects 
prompted  us  to  perform  a  joint  experiment  with  these  colleagues.  The  experiment 
consisted  of  measuring  the  reflection  coefficient  of  a  thermoacoustic  element  as  a 
function  of  frequency.  One  end  of  the  element  was  maintained  at  room  temperature,  the 
other  near  liquid  nitrogen  temperatures.  The  results  of  the  experiment  are  currently  being 
analyzed.19 
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